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ABSTRACT 
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are reviewed. A brief introduction to the status of the experiment is given. The 
review focuses on three areas: multi-scale problems in QCD including results on 
jets, charm, the dynamics of F 2 at low Q 2 , forward jets and fragmentation func¬ 
tions; diffractive structure including results on T production, inclusive diffraction 
in DIS, event shapes and leading baryon production; and, the measurement and 
interpretation of the high-Q 2 neutral and charged current cross-sections. 
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HERA Progress: Since last year’s DIS workshop, the data available for analysis at ZEUS 
has more than doubled, corresponding to an integrated luminosity from the 1994-97 data- 
taking periods of 46.6 pb _1 . This large data sample presents a significant challenge to the 
experimentalists in providing precise data in the many areas which HERA can uniquely access. 
The understanding of the ZEUS detector is being improved using the large data samples to 
calibrate in situ. To illustrate the current level of understanding of the detector: the integrated 
luminosity is currently known to better than ±1.5%; the overall momentum scale of the central 
tracking detector has been established at below the 0.3% level following calibration using elastic 
J/rj; data; the electron energy scale in the barrel calorimeter has been calibrated to ±1% using 
DIS data; and, the hadronic energy scale has been determined within ±3% using DIS data and 
cross-checked using dijet events. Developments in the calibration of the detector, combined 
with the improved statistics, enable increasingly precise as well as new measurements to be 
performed, as discussed below. 


1 Multi-Scale Problems in QCD 


At HERA there is always at least one soft QCD scale, parameterised in terms of the parton den¬ 
sities of the proton which require input from data. Quantities are measured with an additional 
hard scale, such as the photon virtuality, Q 2 , in inclusive F 2 measurements. The measurements 
are then compared to perturbative calculations where the cross-sections may be factorised and 
the soft input constrained. In general, an additional soft (hadronisation) scale enters wherever 
part of the hadronic final state is measured. In addition, in resolved photoproduction processes, 
the photon may also be described as a source of partons introducing an additional soft scale 
which again requires data input. An important point in the development of our understanding 
of QCD is to ensure that the measurements are well-defined theoretically at the required level 
of precision. These multi-scale measurements can then be confronted with the QCD predictions 
to test whether our current understanding is sufficient to describe the process under consider¬ 
ation. A wide range of measurements can be made at HERA which test our understanding of 
QCD and its factorisation properties. 


Photon Structured Before discussing the latest preliminary measurements, it is worthwhile to 
note the uncertainties, as discussed in recently published results on inclusive jet production 
using the cone algorithm.! Experimentally the uncertainty on the jet energy scale of ±3% 
dominates the uncertainty on the cross-sections. Theoretically, the scale dependence on the 
renormalisation and factorisation scales varied between E ^ 1 /4 < fip = Rr < E ^, is minimised 
for a cone radius R ~ 0.6, as observed earlier at the Tevatron. The variation of R sep , the two-jet 
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merging parameter, leads to similar uncertainties on the cross-section for R < R sep < 2 R. In 
order to minimise the theoretical uncertainties due to merging/seed-finding ambiguities, the 
iterative ^-algorithm has been adopted. An inclusive dijet analysis was performed requiring 
E^ 1 > 14 GeV and E^ t2 >11 GeV. At these E^ values, comparisons with Monte Carlo 
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Fig. 1. Inclusive dijet photoproduction cross-sections. 


indicate that multiple interactions (not included in the NLO calculations) are not required 
to describe the observed x 1 distributions or the jet profiles. To gain greater sensitivity to 
the internal structure of the photon the measurements are made at high-y (0.5 < y < 0.85), 
corresponding to the highest accessible photon energies. In Fig. |l] the da/dr} 2 distribution for 
the second-highest energy jet is presented for fixed intervals of r/i. The shaded band represents 
the hadronic energy scale uncertainty. The hadronisation uncertainties have been estimated to 
be ~ 10%, the scale uncertainties are estimated at ~ 10% while the proton parton densities 
are well constrained in the probed x region (x ~ 1CU 2 ). The data are thus sensitive to the 
choice of photon parton densities, as illustrated by the comparison of the the data with the 
GS96 and GRV parameterisations. The ZEUS data has now reached a level of precision where 
the photon parton densities can be discriminated: a global analysis of photon parton densities 
incorporating such data is therefore required. 
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Multijet Structures In order to probe the QCD matrix elements at a deeper level, an inclusive 
three-jet analysis has been performed using the /^-algorithm for jets with > 6 GeV. The 
cross-section for such processes can be written as: 

da( yp jjj) = Yj[J dx^dxp f lh (s 7 ) (x p ) j M 12-^345 1 2 dPS 

M jj i cos 03 , ip 3 X :i ,X 4 
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Fig. 2. Angular distributions for inclusive three-jet events compared to Monte Carlo simula¬ 
tions denoted by the histograms and the 0(aa |) calculations discussed in the text. 


Here, the measured invariant mass distribution, is governed by the photon and proton 
parton densities, /i/ 7 (x 7 ) and fy p (x p ), whereas the scaled energies of the jets, A" 3 and X 4 , are 
controlled largely by phase space. The measured angular distributions in the three-jet centre 
of mass are cos0 3 , where 0 3 is the angle the highest-energy jet makes w.r.t. the beam axis, and 
0 3 , the angle of the three-jet plane w.r.t. the beam axis. The cos0 3 distribution is determined 
by the spin of the primary exchange and the distribution of "03 is related to the coherence 
property of the radiated (lowest-energy) jet to lie in the plane of the beam and the highest- 
energy jet. These are shown in Fig. |2|, compared to 0(aa |) calculations (thick line) as well as to 
a phase space calculation (thin line) where the spin of the partons is ignored. The comparisons 
constitute a refined test of the photoproduction QCD matrix elements Xi 12 ^ 345 - 

Charm in PhotoproductionQ Heavy flavour production introduces a new scale with which to 
test perturbative QCD. Measurements of the D* —» ( K7t)tt s and D* —> (Kh'K'k)tt s channels 
are shown to be in good agreement in Fig. |3](a). Two approaches have been taken in the 
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(a) D* p_i cross-section. 


(b) da/dr) D for p± > 3 GeV. 


Fig. 3. D* ep cross-sections compared to NLO calculations. 


calculations: the “massive” approach (dashed lines) where charm is generated dynamically and 
divergences are regulated by the charm mass; and the “massless” approach (full lines) where 
charm production is activated at the m c threshold and massless approximations are used. 
Uncertainties arise due to the choice of the effective charm mass, and renormalisation and 
factorisation scales (here fiR = \jm 2 c + p 2 ± = /xp/2) as well as the hardness of the charm decay 
to D*, characterised in terms of e c in the Peterson fragmentation function for the “massive” 
calculations. The measured cross-section is typically underestimated in the calculations. In 
Fig. H|(b) this excess is observed to be predominantly in the forward direction. The open 
question is whether the data are more sensitive to the choice of input photon structure function 
or to the D* fragmentation dynamics as we approach the proton fragmentation region or, 
perhaps, anomalously large contributions from b —*■ c decays. Further information is provided 
by the measurement of associated dijets with E^ tx > 7 GeV and Ej7 t2 > 6 GeV shown in 
Fig. [|. Here the measurement of x° BS has a relatively large contribution at low-ay whereas 
the “massive” NLO calculation is significantly more peaked towards one, underestimating the 
low-a; 7 part of the cross-section. The question here is whether the dynamical generation of 
charm in photoproduction is sufficient, assuming that jet hadronisation effects, estimated at 
~ 10% using Monte Carlo simulations, are relatively less important. 

Charm in DIS . i D* —> {K 71)71 s measurements in DIS provide a significant test of the gluon den- 
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Fig. 4. ZU+dijet cross-section versus x° BS 

sity of the proton determined from the scaling violations of F 2 . They will also help to constrain 
theoretical uncertainties in the fits to F 2 where different prescriptions for heavy flavour effects 
are adopted. Compared to the photoproduction case, they remove the uncertainty due to the 
choice of photon PDFs and hence reduce the number of open question posed above. The prelim¬ 
inary cross-section a ep ^ D * x = 8.55±0.40lo;24 nb is measured in the range 1 < Q 1 2 < 600 GeV 2 , 
0.02 < y < 0.7,1.5 < pt(D*) < 15 GeV, and r](D*) < 1.5. In Fig. [|, the upper plots show 
the measurements of the hadronic final state variables pt(D*), r](D*) and Xp*, the fractional 
momentum the D* in the rest frame: the data agree with the massive NLO calculations 
where e c = 0.035, except perhaps at lower xd* corresponding to higher rj(D*). In addition the 
kinematic variables, W, Q 2 and x shown in the lower plots are in good agreement with the 
NLO calculations: it is therefore reasonable to extrapolate the measured cross-section to the 
full {ij(D*)ip T (D*)} range Qto determine F 2 C (±, Q 2 ) via the expression 

In Fig. |6|, the Fg(x, Q 2 ) data mirror the rise of F 2 at small x. The data are in agreement with 
the GRV94 PDFs, where the band represents an estimated theoretical uncertainty due to the 
effective charm mass (m c = 1.5 ±0.2 GeV). This comparison verifies the steep rise of the gluon 
density at low x with a precision of ~ 15 — 20%. 


1 This procedure neglects the possibility of additional contributions outside the measured region due, for example, 

to intrinsic charm. 
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Fig. 5. DIS D* cross-sections compared to NLO calculations (shaded band). 
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Fig. 6. F|(x, Q 2 ) as function of x for fixed Q 2 compared to NLO calculations. 
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(a) A e ff versus Q 2 . (b) dF^/dlogQ 2 versus x. 

Fig. 7. Fits to F 2 data (•) exploring the transition region. 

Transition Region .H The corresponding rise of F 2 with decreasing x, or equivalently the rise 
of c±7 p with increasing W, has stimulated significant theoretical developments in the under¬ 
standing of QCD at high energies. One challenge is to explore how and where the transition 
occurs from soft to hard physics and interpret low-Q 2 data. In order to relate the low-Q 2 
and Q 2 = 0 data, a GVMD (Generalised Vector Meson Dominance) analysis has been per¬ 
formed. This analysis relates the virtual photon cross-section to the real cross-section via 
op! = cbp • A7 ( j /(Ml + Q 2 ), for fixed W (cr L contributions at small Q 2 lead to a small correction 
of a tot ). A good description of the ZEUS BPC data measured in the range 0.1 < Q 2 < 0.65 GeV 2 
is found with Mq = 0.53 ±0.04 ±0.10 GeV 2 . Extrapolating to Q 2 = 0 GeV 2 , the corresponding 
UF 2 ("jp(°) _1 ) dependence is given by the pomeron intercept value 

o±(0)rpc = 1 .145 ± 0.02(stat) ± 0.04(s?/s) (preliminary) 
to be compared with the Donnachie-Landshoff value = 1.08. In this Q 2 range, the rise 

of the cross-section is therefore relatively modest. This behaviour is also seen in the lower Q 2 
points of Fig. [7](a). Here additional datsets are incorporated in fits to the F 2 data of the form 
F 2 = c ■ x XeS | q 2 . The parameter A e fr — a^(0) — 1 for x < 0.01 is then plotted as a function 
of Q 2 . A relatively slow transition from A e fr — 0.1 is observed with increasing Q 2 . This rise of 
F ‘2 with decreasing x is intimately coupled to the scaling violations via the gluon density (in 
leading order dF 2 /dlogQ 2 ~ xg(x) neglecting sea quark contributions). In Fig. 0(b), fits of the 
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form F 2 = a+b-log(Q 2 )\ x have been performed to the published HERA data and the pre-HERA 
prediction from the GRV94 PDFs.@ For x ^ 10 -4 , corresponding to < Q 2 > $ 2 GeV 2 there is 
a qualitative change in behaviour where the scaling violations stabilise and then decrease for 
lower-x values, a behaviour which is not reproduced by the GRV94 PDFs. 

The question is whether this scaling violation behaviour and 
the slow onset of the rise of F 2 with decreasing x can be si¬ 
multaneously understood. A DGLAP NLO fit to the Q 2 > 

1 GeV 2 data (not shown) describes the data, demonstrat¬ 
ing that there is sufficient flexibility in such an approach 
to go down to relatively low Q 2 . However, the relatively 
stable scaling violations observed around <Q 2 >~ 2 GeV 2 
in Fig. [7](b) yield a gluon contribution which is rapidly di¬ 
minishing at small-x and which is significantly smaller than 
the sea quark contribution for small starting scales, as il¬ 
lustrated in Fig. |j: in this low Q 2 region the sea appears 
to be driving the gluon at low-x. For larger Q 2 values the 
gluon dominates the sea and we have an intuitively appealing picture where gluons radiate sea 
quarks. Whether such low-Q 2 partons are universally valid could be tested using e.g. low-Q 2 
data. 

Forward Jet Production .! Why does F 2 rise? In the DGLAP approach the x dependence is 
an input determined at a starting scale Q 2 and evolved in Q 2 . In the BFKL approach the 
x-dependence has recently been calculated in NLO. The underlying dynamics may be tested 
using semi-inclusive forward jet measurements in low-x events. Jets with E\ ~ Q 2 and x jet Jx, 
where Xj et is the momentum fraction of the jet relative to the incoming proton, are selected in 
order to enhance BFKL-like contributions where forward gluons may be emitted at relatively 
large Ft- In Fig. || the jets observed at detector level are shown as a function of E\/Q 2 
compared to three Monte Carlo simulations: LEPTO 6.5 and HERWIG 5.9 are DGLAP-based 
models such that gluons emitted at successively larger x g i uon ~ Xj et have successively lower 
Et whereas the colour dipole model ARIADNE 4.08 incorporates a BFKL feature that gluons 
are not strongly ordered. Three regions are identified: I - E 2 i < Q 2 / 2, the ’’DGLAP” region, 
where all models approximately describe the data; II - Q 2 /2 < E"f < 2 Q 2 , the ’’BFKL” region 
where the DGLAP-based models (LEPTO 6.5 and HERWIG 5.9) fall below the data; and, 
III - E'j, > 2 Q 2 , where all models fail and we may need to describe the DIS data in terms of 
a virtual photon whose structure is being resolved by the high-FV jets. The cross-section is 
evaluated in region II as a function of x for Xj e t > 0.036 and Et > 5 GeV in Fig. [10]. BFKL 
dynamics leads to an enhancement of the forward jet production cross-section proportional to 
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Fig. 8. Low Q 2 parton densities. 
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Fig. 9. Observed E\jQ 2 distribution. 


(. Xj et /x) a p 1 over the 0(aa |) calculation. As shown in Fig. [T0](a) there is a significant difference 
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Fig. 10. Forward jet production cross-sections at (a) parton level and (b) hadron level. The 
data are compared to various predictions discussed in the text. 

between the 0(aa |) MEPJET calculation (represented by the shaded band which includes the 
uncertainty on the renormalisation scale) compared to the leading-order BFKL prediction (full 
curve) at parton level. There are residual uncertainties in determining the hadron-to-parton 
level corrections and therefore the measurement in Fig. p~0](b) is presented at hadron level. The 
rise of F 2 at small x is mirrored by the rise of the measured forward jet cross-section which is 
not described by the DGLAP-based models. A consistent description of the F 2 and forward-jet 
data represents a considerable challenge to our understanding of QCD. 
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Fragmentation Functions: Recently published results on jet shapes have shown that the ob¬ 
served patterns of QCD radiation in liigh-Q 2 neutral current and charged current are similar 
to those observed in e + e _ experiments.! At HERA we are able to study these properties as a 
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(a) ln(l /x p ) max versus ln(Q). 


(b) 1/ada/dxp versus Q 2 . 


Fig. 11. Charged hadron fragmentation in the current region of the Breit frame. 

function of Q 2 in a single experiment and hence provide detailed information on quark frag¬ 
mentation properties.0 These properties can be studied in the semi-soft limit by measuring the 
ln(l /x p ) distributions, where x p = 2 p/Q is the scaled momentum of charged hadrons in the 
current region of the Breit frame. The observed Gaussian distributions are then fitted within 
±1 of the mean to yield the ln(l /x p ) max values as a function of Q given in Fig. |TT|(a). The 
precise data are consistent with data from e + e~ experiments establishing the universality of 
fragmentation over a large range of Q. An MLLA+LPHD fit (indicated by the full line) of 
the form \w{l/x p ) max = 1/2 Y + c 2 \f{Y) — c|, where Y = ln<5/2A e // and c 2 = 0.52 for three 
active flavours in the cascading process, provides a reasonable description of the data with 
A e ff ~ 245 MeV. The high-statistics data enables the region of high x p to be studied. In 
Fig. 0(b) the fragmentation function data are presented for different x v intervals as a func¬ 
tion of Q 2 (in various ranges of x). At high x p and higher Q 2 the fragmentation functions 
exhibit negative scaling violations, consistent with a dominant QCD Compton process (c.f. 
high-a; structure function data). The DIS data (full symbols) are observed to be reasonably 
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Fig. 12. < 1 — T z > versus Q. 


consistent with e + e _ data (open symbols), but systematically lower at intermediate x p values. 
Comparisons with NLO calculations where fragmentation functions extracted from e + e _ data 
are implemented can describe the DIS data (see Fig. 2 in0). However care needs to be taken to 
explicitly include strange-quark fragmentation functions which are systematically softer than 
those of the up/down quarks. This is important since the production of strange quarks from 
the proton sea in DIS is suppressed (by a factor ~ 0.2) compared to those from e + e _ annihi¬ 
lation. In order to compare with recent 1/Q 2 power-correction calculations, the fragmentation 
function data were also presented as function of x\\ = 2 p ■ q/Q 2 = 2 pz/Q, where the Z direc¬ 
tion is defined by the virtual-photon proton axis. Another approach to investigate the role of 
such power corrections is to sum over these momenta and measure the corresponding thrust 
distributions Tz = TLpz/'Zp. A series of event shape variables have been measured in the 
current region of the Breit frame where the power corrections, determined from renormalon 
calculations, are expected to be characterised by a universal aT and to fall as 1/Q. In Fig. [12] 
the measurements using charged tracks in the current region of the Breit frame are displayed as 
< 1 — Tz > versus O.0 The characteristic behaviour is in reasonable agreement with published 
HI results given at slightly different x values. Other event shape variables are chosen which are 
relatively insensitive, in varying degrees, to soft gluon emission and collinear parton branchings 
in order to determine whether a universal aT can be applied. As a first step, comparisons with 
NLO calculations illustrate the need for additional power correction terms (see Fig. 2 in0). 
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2 Diffraction 


The study of diffractive phenomena and their relation to the inclusive cross-section at HERA 
may hold the key to understanding how a “pomeron” structure emerges from QCD. What we 
have learned so far is that a single soft pomeron does not describe all diffractive data measured 
at HERA. As the photon virtuality, Q 2 , or the vector meson mass increases the cross-sections 
rise rapidly with increasing W 2 . As we investigate the pomeron more closely, a new type of 
dynamical pomeron may begin to play a role whose structure is being measured in DIS. These 
inclusive diffractive data are consistent with a partonic description of the exchanged object 
which may be described by perturbative QCD. At this workshop, new results were presented 
on (j) electroproductionjl! the f-dependence of p, (j) and J/-0 photoproduction; 0 ij }'/ J/i\) and T 
photoproduction;0 the determination of and the diffractive/total cross-section in DIS;E1 

diffractive dijet photoproduction cross-sections;@ event shapes in DIS using the LPS;0 and, 
leading baryon production in DIS.0 These results will help to clarify our understanding of pro¬ 
cesses which represent a significant fraction, ~ 10%, of the corresponding total cross-sections, 
yet cannot be understood in terms of a simple hadronisation mechanism. 

Vector meson t dependences .S Diffraction is characterised by a steeply-falling dependence of 
the cross-section as a function of t, the momentum transferred at the proton vertex. This 
characteristic fall-off increases approximately linearly with increasing W. This “shrinkage” 
behaviour is built into the Donnachie-Landshoff pomeron 

a>(t) = a(0) + a' ■ t = 1.08 + 0.25 ■ t. 

Measurements of the quasi-elastic (yp —> Vp) and proton-dissociative (yp -> EA) production 
of vector mesons as a function of t have been performed. The W dependence of the cross-section 
at fixed values of t allows a direct determination of the exchanged trajectory. In order to gain 
sufficient range in W and hence constrain a', the ZEUS preliminary results are combined with 
lower-IU elastic measurements and HI published results. In Fig. [L3|(a) the p° trajectory is 
shown to be a strong function of t with 

a(t) p o = (1.097 ± 0.020) + (0.163 ± 0.035) • t (preliminary) 
for the linear fit where the errors are the combined (stat © sys) errors from each experiment. 
Similarly, the (j) trajectory (where non-leading Regge trajectories are suppressed due to the 
Zweig rule) is determined as 

— (1-083 ± 0.010) + (0.180 ± 0.027) ■ t (preliminary). 

These results are in contrast to the J/^ trajectory shown in Fig. p~3|(b) where 
a(t)j/y = (1.175 ± 0.026) + (—0.015 ± 0.065) • t (preliminary). 

The result for a' is consistent with no shrinkage of the cross-section, which is qualitatively 
expected from perturbative QCD which allows only a slow (logarithmic) dependence of the 
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trajectory as a function of t. These results herald an era of study where the gluon-dominated 
QCD structure at the proton vertex is probed using a diffracted heavy flavour probe. 

ZEUS 1995 preliminary 




(a) p° trajectory. 


(b) J/b trajectory. 


Fig. 13. Pomeron trajectories determined from fits to HERA and lower W data. 


Higher Mass Vector Mesonsx - 3 Twenty years after the discovery of the T, the first observation 
of a signal in photoproduction is shown in Fig. [L4](a) as a broad enhancement around 10 GeV in 
the di-muon mass spectrum. The insert in Fig. |h|(a) indicates the J/ijj and -0' resonances: the 
production ratio is measured to be R = 0.16 ± 0.02 ± 0.04, where the largest contribution 
to the systematic uncertainty is on the tjj' —»branching ratio. The result is in agreement 
with QCD calculations, determined by the wavefunction at the origin for the IS and 2S states, 
of ~ 0.17. The observation of the T leads to the measurement of the elastic cross-section 
for the unresolved Y(1S), T(2 S) and Y(3S) states shown in Fig. 0(b)- Here the relative 
rates or muon production from the T states is determined from CDF data and applied as a 
correction to the T(IS') calculation. Comparison with the leading-order QCD calculations, 
where Odiff ~ xg(x) 2 ®<x, indicates that the measured cross-section is above these expectations. 
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(a) M^+fj- spectrum. (b) T cross-section versus W compared to LO cal¬ 

culations. 


Fig. 14. Measurement of the elastic T photoproduction cross-section. 


Diffractive Structure Functions fl A new era for diffraction was opened by the observation of 
large rapidity gap events in DIS and their subsequent analysis in terms of a diffractive cross- 
section. The diffractive contribution is identified as a non-exponentially suppressed contribution 


at small masses, Mx, of the dissociating virtual photon system. In Fig. 15 the ratio of this 
diffractive contribution to the total virtual photon proton cross-section is given as a function of 
W for various Mx intervals. An approximately constant ratio with W indicates a diffractive 
contribution which rises with a similar W dependence. This simple observation is contrary to 
the naive expectation where the diffractive contribution is identified with the forward part of 
the scattering amplitude and would therefore rise twice as quickly as the total cross-section 
as a function of W. The rise of the diffractive cross-section with W can be parameterised in 
terms of a power law, yielding a p (0) = 1.16 ± 0.01 ± 0.02 (preliminary), after integration over 
t with a mean exponential slope, b = 7 GeV~ 2 and assuming a' = 0.25 GeV~ 2 . In order to 
understand the driving mechanism responsible for this rise, the cross-sections at fixed Mx and 
W are plotted in terms of scaling variables. Integrating over xjp, the momentum fraction of 
the pomeron within the proton, leads to the (3 dependence of (= Fff specified at 

Xp = 0.0042) shown in Fig. [TO] where (3 is the momentum fraction of the struck quark within 
the pomeron. An approximately flat dependence on f3 is observed and an approximate scaling 
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in Q 2 emerges from analysis of the data. The decreasing fraction of events in each M x interval 
with increasing Q 2 observed in Fig. 15 can thus be seen as due to integrating over a decreasing 
(5 ~ Q 2 / (Q 2 + Ml) region which is approximately flat in (3. 


Diffractive Event Shapes . O The measurements of the structure function of the pomeron con¬ 
strain various models of diffraction. These models may be discriminated using event shape 
variables which are directly sensitive to the underlying partonic structure. Similarly, the data 
can be directly compared with e + e _ data where the underlying gluon Bremsstrahlung struc¬ 
ture is well known. Tagging a leading proton in the LPS allows a wide range of M x up to 
25 GeV to be explored for xp < 0.03 and Q 2 > 4 GeV 2 . Measurements of the mean thrust 



Fig. 17. ^Thrust> versus M x 


are presented as a function of M x in Fig. [T7]. The LPS measurements are able to discriminate 
amongst various models which assign different partonic structures to the pomeron. The data 
exhibit reasonably similar values compared to e + e“ data for all values of M x suggesting that 
additional gluon contributions from the pomeron are relatively small. They are also consistent, 
but systematically higher than, the HI data obtained using the large rapidity gap method. 
Improved statistics, from existing data, will help to clarify whether the gluon contribution of 
the pomeron determined from the scaling violations of F^ 2 ' 1 0 is consistent with that extracted 
from event shape variables^ and diffractive dijet photoproduction.0 

Leading B ary on Production . El Non-diffractive contributions play a role at higher values of xp. 
Measurements of forward proton production allow us to compare data with fragmentation 
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models as well as models based on reggeon exchange. Preliminary results for the corrected 
rates of proton production in the range 0.60 < Xl = 1 — xp < 0.91 for DIS at low and higher 
Q 2 (0.1 < Q 2 < 0.8 GeV 2 and Q 2 > 4 GeV 2 ) are (13.0 ± 0.5±g;|)% and (12.7 ± 0.3 ± 0.9)%, 
respectively. The measured rates are typically higher than the currently available fragmentation 
models by a factor of 1.5 to 2. Measurements of uncorrected forward neutron production rates 
have been compared to reggeon exchange models: the comparisons indicate that the internal 
structure of the exchanged reggeon falls like 1 — f3 at large [3 (a dependence similar to n (and 
higher) Regge exchanges), in contrast to the approximately flat behaviour with f3 noted earlier 
for diffractive exchange. The rates for neutron production are approximately the same for DIS, 
photoproduction and even proton-gas interactions although the W dependence of these cross- 
sections is significantly different. This suggests that “diquark” fragmentation is a universal 
process which is to a large extent independent of the type of interaction with the incident 
proton. 


3 High-Q 2 Cross-Sections 

The HERA collider provides a unique window to explore ep interactions at the highest energies, 
extending the range of momentum transfer Q 2 by about two orders of magnitude compared to 
fixed-target experiments. As the HERA luminosity increases we explore the region of Q 2 ~ 
10 4 GeV 2 , where electroweak effects play a role. It is in this unexplored kinematic region that 
we are sensitive to deviations from the standard model (SM). At last year’s DIS97 workshop, 
HI and ZEUS reported an excess of events compared to the SM predictions from the neutral 
current (NC) data taken during 1994 to 1996. For the ZEUS analysis the observed rates 
agreed with expectations except for an excess at the highest Q 2 where two outstanding events 
with Q 2 ~ 40,000 GeV 2 were observed from a luminosity of 20.1 pb _1 . In Fig. |H|(a), these 
events clearly stand out but no new NC outstanding events are observed in the 1997 data, 
corresponding to a further 26.5 pb' 1 of data. In Fig. |l^(b), one charged current (CC) event 
is observed at very high Q 2 ~ 30, 000 GeV 2 from the 1994-1996 data as well as two further 
events from 1997 at Q 2 ~ 20,000 GeV 2 . The number of events is higher than expectations 
but is consistent with the standard model. Attention has therefore focussed on measuring the 
cross-sections at the highest accessible Q 2 values. The theoretical uncertainty on the cross- 
sections was determined from a ZEUS QCD fit to the structure function data on proton and 
deuteron targets from SLAC, BCDMS and NMC as well as the neutrino measurements from 
CCFR, taking into account the correlation amongst the systematic errors of each experiment. 
as(Mz) was varied from 0.113 to 0.123 and a 50% uncertainty in the strange quark content 
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(a) Neutral Current Events. 
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Fig. 18. Reconstructed y versus x scatter plots for (a) neutral current and (b) charged current 
events from the 1994-96 (o) and 1997 (•) data. 


was included. In addition, various published PDFs with different models for charm evolution 
were used as well as fits incorporating E706 prompt photon data and CDF jet data. The results 
yielded SM cross-section uncertainties of ~ 6-8% on the NC cross-section and ~ 6-12% on the 
CC cross-section at the highest accessible Q 2 values. These cross-sections therefore represent 
a benchmark for the standard model. The cross-sections, discussed below, are corrected to the 
electroweak Born level and integrated over the complete y range. 

Charged Current Cross-Sections .0 Charged current events are identified by their missing trans¬ 
verse momentum (pr) due to the escaping neutrino. The cross-section is sensitive to the valence 
d-quark distribution in the proton: 


d 2 cr, 


e^p 


g 2 f 


■[u + c + (l-y) 2 (d + s)j. 


dxdQ 2 2 tt (1 + Q 2 /M 2 v y' 

do cc /dQ 2 was measured for Q 2 > 400 GeV 2 using the Jacquet-Blondel method where Q 2 JB = 
p F /1 — Vjbi with an RMS resolution on Q 2 of ~ 25%, reflecting the 35 %j\[E hadronic energy 
resolution. 869 events were selected with backgrounds below 2.5%, diminishing with increasing 
Q 2 . The systematic uncertainties, mainly due to the hadronic energy scale uncertainty of ±3%, 
correspond to ~ 15% uncertainties on the cross-section at lower Q 2 but increase at larger 
Q 2 . In the upper plot of Fig. [19] the cross-section is observed to fall over more than four 
orders of magnitude. The ratio of the data to the SM, adopting the CTEQ4D PDF, is shown 
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in the lower plot of Fig. |l^ where good agreement is observed up to Q 2 of ~ 10,000 GeV 2 . 
Comparison of the the data uncertainties with those from theory (shaded band) indicates that 
the data will help to constrain the d-quark densities at large-x. The cross-section is suppressed 
at lower Q 2 , due to the 1/(1 + Q 2 /M‘^ v ) 2 propagator term: this characteristic dependence on 
Q 2 has been fitted to yield a value for the mass of the exchanged (space-like) lU-boson of 

M w = 78.6t2A(stat.)+l'^(syst.) GeV (preliminary) 
with an additional PDF uncertainty of ±1.5 GeV. 

Photoproduction of lU-bosons decaying semi-leptonically has been investigated by searching 
for events with a high-pr lepton and missing px with 46.6 pb -1 of data. This is interesting in 
the context of the observed excess of high-pr muons with associated missing px observed by 
Hl.§ In the ZEUS analysis four events are observed in the electron channel where 2.22 ± 0.02 
are expected from W production and 1.24 ± 0.35 from various backgrounds. Similarly, zero 
events are observed in the muon channel where 0.46 ± 0.02 are expected from W production 
and 0.84 ±0.23 from other sources. The ZEUS measurements enable 95%CL limits to be set on 
cr(W)(px lss > 20 GeV) of 2.5 pb and 2.0 pb in the electron and muon channels, respectively. 

Neutral Current Cross- Sections^ High-Q 2 neutral current events are easily identified from the 
high-energy scattered positron. The cross-section is particularly sensitive to the valence w-quark 
distribution in the proton: 

^ + (1 “ -(i-C 1 - V 2 ))*F 3 (x,Q% 

Here, F 2 = _F 2 em (l + $z) is the generalised structure function incorporating 7 ± Z terms which 
is sensitive to the sum of the quark distributions ( xq ± xq ) and xF 3 is the parity-violating 
(Z-contribution) term which is sensitive to the difference of the quark distributions ( xq — xq). 
da NC jdQ 2 was measured for Q 2 > 400 GeV 2 using the double-angle method, with an RMS 
(cr) resolution on Q 2 of ~ 5% (~ 2 — 3%), and cross-checked using the electron method. 
Approximately 38,000 events were selected and the systematic uncertainties, resulting from 
various sources, are typically ~ 2 — 3% at low Q 2 increasing to ~ 10% at higher Q 2 . In the 
upper plot of Fig. ^ the cross-section is observed to fall over more than six orders of magnitude. 
The ratio of the data to the SM, adopting the CTEQ4D PDF, is shown in the lower plot of 
Fig. ^ where good agreement is observed up to Q 2 of ~ 30,000 GeV 2 . Comparison of the the 
data uncertainties with those from theory (shaded band) indicates that the data will constrain 
the parton densities of the proton at large-x. 
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A wide range of new interactions would modify the NC cross-sections in a way which can be 
parameterised by an effective four-fermion (eg —> eg) coupling. Given a convention for the 
strength of the coupling (g 2 = 47 t) , limits can be placed on the effective mass scale (A) of these 
contact interactions. Scalar and tensor terms are constrained by earlier experiments and atomic 
parity violation experiments provide strong constraints on various vector couplings: the relative 
size and sign of individual terms in the contact interaction amplitudes is therefore limited to 
24 different combinations. These contact interactions all contain a term proportional to 1/A 4 
which enhances the cross-section as well as a SM interference term proportional to 1/A 2 which 
can either enhance or suppress the cross-section at intermediate Q 2 . No significant deviations 
are found and limits on the 24 models are set in the range of A ~ 2 — 5 TeV. These limits are 
competitive with, and in some cases extend, those limits set from from hadronic cross-section 
measurements at LEP and Drcll-Yan electron pair production at the Tevatron. 


Conclusions 


Selected highlights from the ZEUS analyses of HERA data have been presented. As G. Sterman 
noted at the DIS97 workshop§l “We must use the QCD we know well to investigate new physics; 
but we must also pursue the QCD we do not know well.” At the DIS98 workshop a wealth of 
precise ZEUS data enabled us to identify weak links in our understanding of QCD, to start to 
probe the electroweak sector and to set limits on physics beyond the standard model. 
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